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Abstract 

Hydrogen single-particle dynamics in solid LiH at T = 20 K has been studied through the 
incoherent inelastic neutron scattering technique. A careful analysis of the scattering data has 
allowed for the determination of a reliable hydrogen-projected density of phonon states and, from 
this, of three relevant physical quantities: mean squared displacement, mean kinetic energy, and 
Einstein frequency. In order to interpret these experimental findings, a fully-quantum microscopic 
calculation has been carried out using the variational Monte Carlo method. The agreement achieved 
between neutron scattering data and Monte Carlo estimates is good. In addition, a purely harmonic 
calculation has been also performed via the same Monte Carlo code, but anharmonic effects in H 
dynamics were not found relevant. The possible limitations of the present semi-empirical potentials 
are finally discussed. 

PACS numbers: 61.12.-q,67.80.Cx,63.20.-e,02.70.Ss 
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I. INTRODUCTION 



Hydrogen forms stable stoichiometric hydrides by reaction with all of the alkali metals: 
Li, Na, K, Rb, and Cs.— Since the first x-ray study^ in 1931, it has been shown that LiH, NaH, 
KH, RbH and CsH (AlkH in short) crystallize with the rock-salt structure (i.e. a fee with a 
two-atom basis: Li in (0, 0, 0) and H shifted by (a/2)(i-|-j-|-k)) at room temperature. In these 
materials, experimental evidence seems consistent with hydrogen being present in the form 
of anions or modified anions: electron distribution investigations^ estimated the ionic charge 
in LiH to fall in range from 0.4 to 1.0 electron-charges, indicating that the alkali hydrides 
are probably very similar to the alkali halides, with respect to the electronic structure, and 
that AlkH might even be regarded as the lightest alkali halides, not far from alkaline fluoride 
compounds: AlkF. In addition, calculations of the electronic charge density in LiH, based 
on the local density approximation, suggest that most of the electronic charge is transferred 
from Li to H.- Because of this fact, which gives rise to long-range interactions between 
hydrogen atoms, the H dynamics in these compounds is expected to be very different from 
the one in group HI- VIII metal hydrides. However, apparently only LiH has been studied 
in some detail,— both theoretically and experimentally. 

The choice of lithium hydride (and deuteride) has not been casual at all: they are rock- 
salt crystals having only four electrons per unit cell, which makes them the simplest ionic 
crystals in terms of electronic structure. Moreover, there is a large isotopic effect provided 
by the substitution of the proton with the deuteron. Last but not least, owing to the 
low mass of their constituent atoms, these two compounds represent a good case for the 
calculations of the zero-point motion contribution to the lattice energy. In this respect, 
they are expected to behave as quantum crystals, somehow related to the hydrogen (H2 and 
02)^ or the helium ('^He and ^He)^ families. By comparing an important anharmonicity 
parameter, namely the zero-temperature Lindemann ratio® to the ones of other quantum 
crystals, one finds that H in LiH {'Jh = 0.12)3 lies in between solid H2 (7 = 0.18) and solid 
Ne (7 = 0.09), placed very close to solid D2 (7 = 0.14). Because of these unique physical 
properties (and incidentally because of their use in thermonuclear weapons) LiH and LiD are 
in general fairly well described, even in condition of very high pressure {p > 10 GPa): both 
structurally (neutron diffraction),— and dynamically (second-order Raman spectroscopy),— 
and through ab-initio electronic structure simulations.— 



2 



However, as far as complete (i.e. including full phonon dispersion curves) lattice dynamics 
works are concerned, the situation looks much less exhaustive: LiD dispersion curves have 
been measured by coherent inelastic neutron scattering in 1968,— even though the longitu- 
dinal optical branch was almost unobserved. These experimental data were fitted through 
a 7-parameter shell model (SM7) and converted into equivalent LiH data. Later on,«14 the 
same data have been fitted again through a more advanced method: a deformation dipole 
model with 13 adjustable parameters (DDM13), which also provided elastic and dielectric 
constants, effective charges, second-order Raman spectra, etc., all in good (or fair) agree- 
ment with the known experimental values. The only noticeable exception was represented 
by the H-projected density of phonon states (H-DoPS) for LiH:— 

1 ^ 

^h(w) = — 5mieH(q,j)l^'^(^-^(q,i)) , (1) 

-^^V q j = l 

where q is a phonon wave-vector contained in the first Brillouin zone, N is the number 
of these wave-vectors, j is labeling the six phonon branches, eH(q, j) is the polarization 
vector for H, and uj{q^,j) is the phonon frequency. LiH H-DoPS turned out both from 
SM7 and DDM13 to be rather different from the old incoherent inelastic neutron scattering 
(IINS) measurement In addition, both SM7 and DDM13 are force constant models and do 
not provide any suitable inter-ionic potential scheme for LiH. The potential approach was 
first attempted by Hussain and Sangster,— who tried to include alkali hydrides in a larger 
potential scheme derived for all the alkali halides making use of the Born-Mayer functional 
form: considering the little number (i.e. 3) of adjustable parameters, the LiD dispersion 
curves obtained from this potential are quite in good agreement with the coherent neutron 
scattering data.-^^ However, the existence of some problems in this scheme for LiH (e.g. the 
H-H Pauling parameter turned out to be unphysical, the elastic constants were too large, 
etc.), prompted Haque and IslamA^ to devise a new set of potentials (HI) for both LiH and 
NaH. But, despite their better description of some of the macroscopic properties of LiH and 
their superior physical soundness, the HI potential does not provide any advantage over the 
older one as far as the dispersion curves are concerned. Finally, it is worth mentioning an 
important ab-initio calculation of the LiD dispersion curves,^^ which correctly included the 
zero-point effects in the lattice energy minimization procedure. 

As for thermodynamics, accurate constant-pressure specific heat measurements on LiH 
and LiD in the temperature range T = (4 — 300) K are reported by Yates and co-workers.— 
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These experimental results, once corrected for thermal expansion with the help of tempera- 
ture dependent lattice parameters,— were used to extract the effective Debye temperature, 
Qd, as a function of T in the aforementioned temperature range. The relative thermal 
variation of O/j in LiH was found to be very large: 20% between 4 K and 80 K; more than 
30% between 80 K and room temperature, revealing a strong anharmonicity in the LiH 
lattice dynamics. Later, a comparison was made between experimental Qd{T) data and 
the DDM13 results of Dyck and Jex.-'''* The agreement was found quite good between 30 K 
and 300 K, while for temperature values below 20 K some anomalies in the behavior of the 
experimental specific heat were found, i.e., the measured heat capacity displays a peak at 
Tc = 11.1 K (for LiH), the origin of which has not yet been discovered. The possibility of a 
phase transition to a CsCl-type structure has been proposed,— but it seems rather unlikely. 
The specific heat measurements in the reported temperature range are mainly sensitive to 
the details of the phonon density of states in its low-energy region (say — 30 meV),— which 
lies for LiH in the acoustic band, spanning the phonon energies in the (0 — 70) meV range. 
Thus, unfortunately, thermodynamics seems rather unable to probe the real H dynamics in 
LiH. 

Given the aforementioned scenario, a study on the hydrogen vibrational dynamics in LiH 
has to answer a number of questions that we have tried to synthesize in four main points: 
1) How is the H vibrational dynamics in LiH described by the proposed force-constant and 
potential schemes? 2) How do they compare with new and reliable incoherent inelastic 
neutron data? 3) Are low-temperature anharmonic effects (i.e., purely quantum-crystal 
effects) detectable in the H dynamics in LiH? 4) In case, how is it possible to calculate them 
using microscopic theory? We believe that a combined use of incoherent inelastic neutron 
scattering from which the H-DoPS can be worked out, and of fully quantum simulations, from 
which important equilibrium physical quantities can be derived (mean squared displacement, 
mean kinetic energy and Einstein frequency) can provide a new and deep insight into the 
problem of quantum hydrogen dynamics in condensed matter. The rest of the present work 
will be developed according to the following scheme: the neutron measurements will be 
described in Sect. II, and it will be also shown how to extract a reliable H-DoPS from the 
experimental spectra. Sect. HI will be devoted to describe the simulation technique, namely 
the variational Monte Carlo method. Then, in Sect. IV the results on the H dynamics in LiH 
will be discussed. Here a comparison between the quantities derived from the experimental 
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spectra and their estimates obtained through the variational Monte Carlo simulations, and 
other less advanced techniques, will be established. Finally, in Sect. V we will draw some 
conclusions. 

II. EXPERIMENTAL AND DATA ANALYSIS 

The present neutron scattering experiment was carried out on TOSCA-II, a crystal- 
analyzer inverse-geometry spectrometei^ operating at the ISIS pulsed neutron source 
(Rutherford Appleton Laboratory, Chilton, Didcot, UK). The incident neutron beam 
spanned a broad energy {E) range and the energy selection was carried out on the secondary 
neutron flight-path using the (002) Bragg reflection of 7 graphite single-crystals placed in 
back-scattering around 137.7°. This arrangement fixed the nominal scattered neutron en- 
ergy to E' = 3.32 meV. Higher-order Bragg reflections were filtered out by 120 mm-thick 
beryllium blocks cooled down to a temperature lower than 30 K. This geometry allows to 
cover an extended energy transfer (3 meV< huj < 500 meV) range, even though the fixed 
position of the crystal analyzers and the small value of the final neutron energy imply a 
simultaneous variation in the wave-vector transfer, Q. In other words, on TOSCA-II Q is 
a monotonic function of lu: 2.8 A^^ < Q{uj) ~ cj^/^ < 16.5 A^^. The resolving power of 
TOSCA-II is rather good in the accessible energy transfer range: Ahu/E ~ 1.3 — 2.3%. 

The sample cell was made of aluminum (0.8 mm-thick walls) with a squared-slab ge- 
ometry, exhibiting an internal gap of 1.2 mm. The cell area (85.0 mm x 85.0 mm) was 
rather larger than the neutron beam cross-section (squared, roughly 40 mm x 40 mm). The 
scattering sample was made of 5.733 g of polycrystalline LiH (powder from Sigma- Aldrich, 
97% assay) and contained purely natural lithium (92.5% '^Li and 7.5% ^Li). After collect- 
ing background data on the empty cryostat, we cooled down the empty sample cell to the 
experimental temperature (T=20 K), starting to record the time-of- flight (TOF) neutron 
spectrum as T < 30 K, up to an integrated proton current of 344.1 fiA h. Then the system 
(LiH sample + can) was loaded into the cryostat, cooled down to 20.0 K and measured up 
to an integrated proton current of 2101.0 fiA h. The thermal stability of this measurement 
was quite satisfactory, since the temperature fluctuations never exceeded 0.3 K, and the 
temperature gradient between top and bottom of the cell was lower than 0.1 K. The mean 
sample temperature was estimated to be (20.1 ±0.1) K. Special care was devoted to prevent 
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possible lithium hydroxide formation during the sample loading procedure. In addition, 
the subsequent inspection of the raw spectroscopic data exhibited no sharp features in the 
(75-78) meV range, which are a clear signature of lithium hydroxide contamination.— 

The raw neutron spectrum exhibited a series of strong features (in the 73 — 140 meV 
range) related to the density of optical-phonon states and essentially due to the H and Li 
anti-phase motion in the lattice unit cell. The first overtones of these optical bands are 
also well visible in the recorded spectrum at about twice the aforementioned energy transfer 
intervals. On the contrary, at low energy transfer huj, the acoustic band appears rather 
weaker than the optical one, since here lithium and hydrogen ions move basically in phase, 
and the H mean square displacement is much smaller, more than 20 times according to 
lattice dynamics simulations^^. 

The experimental TOF spectra were transformed into energy transfer data, detector by 
detector, making use of the standard TOSCA-II routines available on the spectrometer, and 
then added together in a single block. This procedure was justified by the narrow angular 
range spanned by the detectors, since the corresponding full-width-at-half-maximum, A^, 
was estimated to be only 8.32° (see Ref. 3])- In this way, we produced a double-differential 
cross-section measurement along the TOSCA-II kinematic path {Q{uj),uj) of the (LiII+ can) 
system, plus, of course, background and empty can. Then, data were corrected for the k' / k 
factor, and background and empty-can contributions were properly subtracted. At this 
stage the important corrections for the self- ah sorption attenuation and multiple scattering 
contamination were performed, the former being particularly important due to the large 
natural Li absorption cross-section (craf,(Li) = 70.5 barn).— These two corrections were 
applied to the experimental data through the analytical approach suggested by Agrawal 
in the case of a fiat slab-like sample. We made use of the simulated H- and Li-projected 
densities of phonon states derived by Dyck and Jex— in order to evaluate: 1) the H total 
scattering cross-section, known to be largely dependent on cj; 2) the LiH scattering law to be 
folded on itself in order to generate the multiple scattering contributions. Both procedures 
were accomplished in the framework of the incoherent approximation,^^ totally justified by 
the preponderance of scattering from H ions, and by the polycrystalline nature of the sample. 
Multiple scattering was found to be around 6.7% of the total scattering in the energy transfer 
range of main interest (i.e. optical phonon region): 50 meV< hu < 150 meV, and then 
subtracted as shown in Fig. Q (a)- After performing the two aforementioned corrections. 
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our neutron spectrum still contained some scattering from the Li ions, which was however 
much lower than the one from the H ions, and, moreover, mainly localized in the acoustic 
phonon region , ^'^i^^ i.e. for huj < 50 meV. So the Li contribution was simulated (always in 
the framework of the incoherent approximation) through the Dyck and Jex calculations of 
the Li-projected densities of phonon statesi^ and then removed. All the practical details 
of this procedure can be found in Ref. 23], where this is applied to various binary solid 



systems, such as H2S, D2S, and HCl, measured on TOSCA-I. 

The last stage before the extraction of the H-DoPS was the evaluation and the subtraction 
of the multiphonon contribution, not totally negligible because of the Q-values attained by 
TOSCA in the 3 meV < hcu < 150 meV range (namely 2.8 A"^ < Q < 9.5 A"^). Processed 
LiH data, proportional to the self inelastic structure factor^^ for the H ions, Ss{Q,uj), were 
analyzed through an iterative self-consistent procedure^^^ aiming to extract the one-phonon 
component of Ss{Q,uj) (see Fig. [T] (b)) and the hydrogen Debye- Waller factor, namely 
S^{Q,u!) and W{Q), respectively. Once again, all the technicalities can be found in Ref. 
j22|, while in the following only the final equation of the procedure is reported. From 
Ss,i{Q,uj) and W{Q), given the isotropic nature of the LiH lattice, Zn^u) is simply worked 



out via!^ 



coth — + 1 

\2kBT 



(2) 



hQ 

where mn is the proton mass. The result for Zi{{uj) is plotted in Fig. |21 Equation and 
the other formulas^l used to work out Zh(co'), are formally exact only in the framework of 
the harmonic approximation, but they have a practical validity that is far more general, as 
explained by Clyde'' in the context of solid helium. From Zh(co'), making use of normal and 
Bose-corrected moment sum rules,— we are able to derive three important quantities related 
to the hydrogen dynamics in this hydride, namely the H mean squared displacement (u^), 
the H mean kinetic energy (Th), and the H Einstein frequency, Qo,}i'- 

(u^) = / dco ^ ' coth ' * 



2mii Jo uo \2kBT J 

^OH = / dujZY{{uo) Up . (3) 

JO 

As for the first, we found: (u|) = 0.062(1) A^, while for the second: (Th) = 80(1) meV, 
and finally for the last: /il^o.H = 109.2(9) meV. 



III. VARIATIONAL MONTE CARLO SIMULATION 



We have studied the ground-state properties of sohd LiH by means of the variational 
Monte Carlo (VMC) method.— VMC is a fully quantum approach which relies on the vari- 
ational principle; it has been extensively used in the past in microscopic calculations of 
quantum crystals, mainly ^He and ^He. Assuming only pair-wise interactions between the 
different atoms in the crystal, the Hamiltonian describing LiH is: 

where Vij is the distance between the atoms composing an i,j pair; mu is the Li average 
atomic mass; A'h and Nu stand for the number of H and Li atoms, respectively; and V^^'^\ 
\/(Li,H) y{u,u) Yepiesent the three pair-wise interaction potentials. The variational 
principle states that for a given trial wave-function \E', the expected value of H is an upper 
bound to the ground-state energy Eq: 

The multidimensional integral required in the calculation of Eq. (jSj) can not be evalu- 
ated exactly by analytical summation methods like, for example, the hypernetted chain 
formalism.i^ However, a stochastic interpretation of this integral is rather straightforward, 
and this is actually the task carried out by the VMC method, with the only cost of some 
statistical noise. 

A key point in the method is the search for a trial wave-function with a sizeable overlap 
with the true ground-state wave-function \E'o- An extensively tested model in quantum 
crystals is the Nosanow-Jastrow trial wave-function.-^'^ According to this description, the 
wave-function is written as a product of a Jastrow factor F, accounting for the correlations 
induced by the interatomic potentials, and a phase dependent term $, which introduces the 
crystal symmetries in the problem: 

^ = F$ . (6) 

The Jastrow factor contains two-body correlation functions between the different pairs f^^'^\ 

j(Li,Li) g^j^^ j(H,Li). 

i<j i<j i,j 



and $ localizes each particle around the lattice equilibrium sites of the crystal phase Rf 
through the functions g^^^ and g^^^^: 

Nil Nu 

<^ = I[9^''\\r. - Rf I) n^^^'^Hlr. - R.^i) . (8) 

i i 

The VMC simulation is carried out at a molar volume vq = 10.059 cm^, which is the 
experimentally estimated value at zero pressure and zero temperature j^^*^ The lattice con- 
stant, derived from the experimental molar volume Vq, is a = 4.0578 A. The calculation 
is worked out with a simulation cubic box containing 108 particles of each type with peri- 
odic boundary conditions. We have checked that this number of particles is large enough 
for practically eliminating size effects on the more relevant quantities in the present study: 
the H mean kinetic energy, its mean squared displacement around the lattice sites, and its 
Einstein frequency. 

As we have seen in the introductory section, LiH and the rest of hydrides and deuterides 
of light alkali metals are well described as ionic crystals. This nearly perfect ionic bonding 
allows for a model in which ions H~ and Li^, both with a Is^ electronic configuration, are oc- 
cupying the lattice sites of the crystal. The present simulation relies on this model, assuming 
rigid closed-shell ions interacting via central interatomic potentials. The overlap repulsion 
potential between the ions is taken from the semi-empirical Born-Mayer interaction,— and 
van der Waals attractive terms are also included to deal with polarizability effects. Accord- 
ing to this general scheme, as previously mentioned, Haque and Islam^^ proposed their pair 
potential V^''^\r): 

viT'\r) = B^, exp {-A^,r) - ^ , (9) 

with {a, (3} = Li, H. The set of parameters entering V^i'^\r) is reported in Tab. 1. In order 
to evaluate the influence of the interatomic potentials in our results we have also used the 
model proposed by Sangster and Atwood,— Vg^'^''(r): 

^s^'^^ (r) = B^, exp {E^p + F^, - r)] - ^ - ^ , (10) 

which incorporates a dipole-quadrupole attractive term. The parameters for LiH, which are 
taken from Hussain and Sangstei*ii (see also Sect. I), are reported in Tab. 11. 

The main task in a variational approach like the present one is to seek for a good trial 
wave-function \E' (Eqs. (jEHHI))- Like in solid helium calculations,^^ we have chosen analytical 



9 



two-body /(°^'/^)(r) (see Eq. ((7j)) and one-body g^"\r) (see Eq. (jH))) correlation factors with 
a set of free parameters to be optimized. In particular, /(°''5)(r) is of McMillan type?^ 

5" 



/("'^)(r) = exp 



1 (bqp 

2 [ r 



(11) 



and the specific phase factor is a Gaussian centered on the sites of the crystal lattice: 



^?(")(r) = exp(-^c,r2) . (12) 

In the optimization search the absolute minimum of the internal energy {H)<f, is looked 
for. However, the number of parameters is large enough to make this optimization rather 
difficult. In order to discern between local minima of similar quality we have also calculated 
the energy of the solid including Coulomb contributions, considering the ions as point-like 
particles. Under this criterion, the final parameter set is the one which simultaneously 
minimizes the short-range energy and the total energy including Coulomb contributions. 
The values obtained, which are the same for the two short-range potentials Vhi and Vsa, 
are: Cu = 150 A"^ cr = 15 bun = 2.0 A, buu = 1-5 A, and 6hh = 1-0 A. Among 
the three ba/B Jastrow parameters, the most important is the cross one {bun), since the first 
neighbors of a H~ ion are Li"*", and vice versa. It is worth noticing that the optimal values 
for Cq, refiect the difference in the degree of localization around the sites between the two 
ions due to their significantly different masses, ch << cu- Including the Coulomb potential 
in the energy, the energies per particle are -5.07(2) eV and -5.32(2) eV for the Vm and Vsa 
short-range potentials, respectively. 

Information on the spatial structure of the solid can be drawn from the two-body radial 
distribution functions g^°''^\r). In Fig. 3, results for the three components g^°''/^^(^r) are 
shown. As expected, the location of the peaks follows the inter-particle pattern imposed 
by the lattice: each ion is surrounded by ions of opposite sign and ions of equal sign are 
distributed with the same periodicity. The major mobility of II~ with respect to Li+ is also 
observed by comparing the height and the spreading around the sites of g^^^'^^^r), on one 
hand, and (y((^'^)(r) and to a lesser extend g^^^'^\r), on the other. 

A structural quantity which can be directly compared with the present experimental data 
is the mean squared displacement of the H~ ions, (u^): 

>ua = ^(f:(^.-Rr)V (13) 



\i=l 
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Using the trial wave- function quoted above, the VMC result is (u^) = 0.074(2)12. rpj^^ 
resulting Lindemann ratio: 



vK) ,,,, 

7H = , (14) 

is 0.134. Additional insight on the spatial localization of H" can be obtained by calculating 
the H~ density profile mh('"), with r being the distance between the ion and its site. The 
function Mh('") is shown in Fig. 4; it is very well parameterized by a Gaussian with the VMC 
expected value (u^) = 0.074(1) (see the solid line in the same figure). 

The kinetic energy per particle of H~ is one of the partial contributions to the total energy 
of the solid, which is evaluated at each step of the VMC simulation. It is the expected value 
of the operator: 

with configuration points generated according to the probability distribution function 
As in the estimation of (u^), the H~ kinetic energy is the same for the two interatomic 
potentials (Eqs. I9I10|) since the optimization procedure has led to the same variational 
wave-function. The result obtained is (Th) = 84(l)meV. 

The third physical quantity evaluated in the present experiment, the Einstein frequency 
r2o,H; can be calculated in a VMC simulation through its proper definition: 

K.-^^(plV^, (16) 

where the expected value of V^Vh(t') is calculated over the configurations generated by the 
probability distribution function and Vh('") is the potential felt by an H~ ion: 

Vnir.) = \ E ^^"'"Hr.) + ^ E V^^^.^ . (17) 
The result obtained applying Eq. ([1^) is TiVLq yi = 110.3(5) meV. 



IV. DISCUSSION 



The aim of this discussion section is threefold: 1) critically analyzing various dynamical 
quantities derived from the scientific literature with respect to the present results obtained 
through neutron scattering (Subsect. A); 2) comparing the outputs of the experimental and 
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theoretical approaches presently used, namely TINS and VMC simulations (Subsect. B); 3) 
finally, shedding some light on the important point of the possible quantum anharmonic 
effects on the H~ dynamics in low-temperature LiH (Subsect. C). 

A. Analysis of the H-DoPS evaluations 

A comparison among the various determinations (both experimental and simulated) of 
the hydrogen-projected density of phonon states could be finally established at this stage. 
In Fig. 121 four H-DoPS estimates have been plotted together in the frequency region con- 
cerning the two optical bands (65 meV < hu < 160 meV), namely the translational optical 
(TO) and the longitudinal optical (LO), which actually contain more than 97% of the total 
Zh(co') area.— The present IINS experimental result is plotted as circles, the full line repre- 
sents the DDM13 lattice dynamics simulation,— the dashed line is the SM7 lattice dynamics 
simulation,i^ and the dotted line stands for the old IINS measurement.-^ As a preliminary 
comment, one can easily observe the existence of a fair general agreement among all the four 
Zii{uj) in the TO range, at least as far as the peak position is concerned. On the contrary, 
the LO region looks much more uncertain, the peak centroid varying from 115 meV up to 
140 meV. The reason for such a behavior is easily understandable for SM7 and DDM13 data: 
these lattice dynamics calculations made use of the some parameters (7 and 13, respectively) 
derived from a fit of the same LiD dispersion curves measured by Verble^'^ (for a detailed 
comparative discussion on the differences between the SM7 and DDM13 H-DoPS calcula- 
tions see Ref. Q|). By a simple inspection of these experimental dispersion curves, it is 
clear that the LO neutron groups are really few (four values plus one infra-red measurement 
at the r point). However the disagreement between the present IINS Zi{{uj) and the old 
one is difficult to explain, so that we are inclined to think that these discrepancies are due 
to experimental imperfections in the data analysis of the latter (e.g. multiple scattering or 
multiphonon scattering subtraction). Selecting the two most recent experimental and simu- 
lated Zh(co'), namely the present IINS and the DDM13 estimates, we can observe an overall 
semi-quantitative agreement, the main discrepancies being concentrated in two regions: at 
low energy, in the onset of TO band (65 meV < huj < 90 meV), and in the LO band as 
a whole (112 meV < hu < 145 meV). As for the latter, a simple energy shift of 4.5 meV 
seems sufficient to largely reconcile IINS and DDM13, while in the former case, neutron data 



12 



appear somehow broader than lattice dynamical ones (TINS FWHM being about 5.1 meV 
larger than DDM13 FWHM). Considering the TOSCA-II energy resolution in this region 
{Afiuj ^ 1.6 meV), a simple explanation based only on experimental effects can be easily 
discarded. However, as pointed out by Izyumov and Chernoplekov for other hydrides,— 
such a broadening of the H-DoPS TO bands might be the mark of the hydrogen anharmonic 
dynamics in LiH, through a finite phonon life-time. In this respect, more will be said in 
Subsect. C. 

B. Comparison between IINS and VMC results 

Since VMC, as we have seen in Sect. HI, is a ground-state method, Zuiuo) can not 
be directly evaluated. However, through the aforementioned normal and Bose-corrected 
moment sum rules in Eqs. 0, it was possible to describe the main features of the H-DoPS 
via (u^), (Th) and ^o^h, which are equihbrium quantities calculated by the VMC code 
(see also Sect. HI). Before proceeding with this comparison, it is worth noticing that the 
VMC calculation is performed at zero temperature, whereas the measure is accomplished at 
T = 20.1(1) K. However, thermal effects are negligible since the Debye temperature of LiH 
is approximately 1100 K,^ and therefore the measured system can be certainly considered in 
its ground state, at least as far as the H~ ion dynamics is concerned. This assumption can 
be easily proved by calculating (always from the experimental Zi{{uo)) the zero-point values 
of the H mean squared displacement and mean kinetic energy, setting T = 0: (u^)(T = 
0) = 0.062(1)1^ and (Th)(T = 0) =80(1) meV, identical within the errors to the values 
estimated at T = 20.1 K in Sect. II. 

Going back to VMC, one can notice a value of the zero-point H mean squared displacement 
(0.074(2) A^, as in Sect. HI) slightly higher than the IINS experimental measure. In addition 
these two figure have to be compared to the most recent neutron diffraction estimate by 
Vidal and Vidal-Valat>a (u|) =0.0557(6) (extrapolated at T = 20 K by the present 
authors from the original data in the temperature range 93 K— 293 K), which appears close 
but still discrepant from the VMC and IINS findings. However, it has to be pointed out 
that a previous room-temperature diffraction result by Calder et al.,'^ (u^) =0.068(1) A^, 
seems to exhibit a similar trend if compared to the Vidal and Vidal- Valat's figure in the 
same conditions: (u^) =0.0650(6) A^. On the other hand, for the other two aforementioned 
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physical quantities, the agreement between IINS ((Th) = 80(1) meV, ^Oo,h = 109.2(9) meV) 
and VMC ((Th) = 84(1) meV, /?.fio,H = 110.3(5) meV) is much more satisfactory, confirming 
the vahdity of our combined IINS- VMC method. 

An interesting test on the obtained results can be accomplished in the framework of an ap- 
proximate estimation known as the Self Consistent Average Phonon (SCAP) formalism .^^i^^ 
The SCAP approach relies on the well-known Self Consistent Phonon (SCP) method, but 
replacing the sums of functions of the phonon frequencies by appropriate functions of an 
average-phonon frequency. Results^ for quasi-harmonic and harmonic solids like Ne, Kr, 
and Xe obtained using SCAP have shown good agreement with experimental data. The 
application of this formalism to quantum crystals seems however more uncertain due to the 
relevant increase of anharmonicity. As LiH seems to be a quantum crystal, but with less 
anharmonicity than for example ^He, SCAP can somehow help in the present study. Nor- 
mally SCAP is used to evaluate various physical quantities in an iterative way, employing 
lattice parameters and interatomic potentials only.— Here, on the contrary, the method will 
be applied in one single step, starting from "exact" values of I^ch- To this end, we have 
calculated via SCAP (at T = 0) the H mean squared displacement and mean kinetic energy 
using the relations: 

3h 



2 \(SCAP) 



(Th)(^^^^) = I mu{ul)^'^^^^nln = ^^^o,h • (18) 

The results obtained through this approximation are: (^^^^{scap-vmc) _ 0.0564(3) A^, 
(u2j)(scAP-iiNS) ^ 0.0570(5) A2, (Th)(scap-vmc) ^ 82.7(4) j^eV, and {T^)iscAP-ims) ^ 

81.9(7) meV. By comparing these approximated values with the microscopic ones, one re- 
alizes that the H mean kinetic energy values come out very close, but the SCAP values 
of the mean squared displacement are significantly smaller, actually not far from the neu- 
tron diffraction result by Vidal and Vidal-Valat.^ The physical meaning of these results is 
straightforward: both SCAP equations are exact at T = in presence of a purely harmonic 
Einstein solid (i.e. if Zh(co') =6 {uj — Qq^u) )• But if the solid system exhibits a broader 
H-DoPS, then (u^)*^^*-^^^) comes out rather underestimated, since f^Qjj is exactly computed 
stressing the high-frequency part of Zh(u;) via the integrand factor (as in Eqs. ©). 
In this way, one somehow corrects for this bias by expressing the zero-point mean kinetic 
energy as the product of (uh)*^^*"^^^ times f^o.H; because (Th) would be exactly evaluated via 
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the integrand factor oj (see Eqs. (jS))), which still stress the high-frequency part of Z'h(co'), 
but less than ^q^. Nevertheless, it is quite remarkable the good accuracy achieved by this 
relatively simple approach in evaluating the mean kinetic energy, probably beyond what is 
a priori expected for a system with a possible quantum- crystal character like LiH. 

C. Possible quantum anharmonic efFecs in LiH 

Given the relatively large value of the H Lindemann ratio in LiH (7// ~ (0.12 — 0.13) 
as seen above), it is natural to inquire on the low-temperature anharmonic effects in the 
ion dynamics. In this respect our tools are well suited, since VMC is a microscopic 
quantum simulation, assuming neither the harmonic approximation like the usual lattice 
dynamic calculations, nor the semi-classical treatment of the particle motion like the molec- 
ular dynamics approach.'^ The method applied to test the existence of possible quantum 
anharmonic effects in H dynamics (at T = 0) was simply devised running the same VMC 
code in a "harmonic way", i.e., replacing for every atomic pair z, j the exact pair potential 
value V^("'^)(ri^) hjM 

where and Uj are the instantaneous position of an i atom and its displacement from 
the equilibrium position, respectively, while r^j stands for the vector rj — r^, and r^^ij is the 
equilibrium separation of an atomic pair i, j. It is worth noting that both the static potential 
energy V^*^"'^''(ro,jj) and the Hessian components are all calculated only once per each atomic 
pair z, j, since they depend only on the equilibrium distance r^ij. However the "harmonic" 
results from VMC did not show any significant difference (within their uncertainties) from 
the "exact" ones (see above in Subsect. B), proving that, at least in the framework of the 
semi-empirical pair potentials employed, quantum anharmonic effects are totally negligible 
in the evaluation of (uy) and (Th), even at T = 0. 

V. CONCLUSIONS AND PERSPECTIVES 

Dynamical properties of solid LiH at low temperature have been studied using incoherent 
inelastic neutron scattering with higher accuracy than in previous measurements. The anal- 
ysis of the scattering data has allowed for the extraction of a reliable hydrogen-projected 
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density of phonon states. From this physical quantity it was possible to estimate three 
relevant quantities intimately related to the microscopic dynamics of H in the solid: its 
mean squared displacement, mean kinetic energy, and Einstein frequency. Apart from the 
intrinsic interest in an accurate quantitative determination of these magnitudes, we have 
tried to shed some light on two fundamental questions on the physical nature of LiH, i.e. its 
quantum character and the degree of anharmonicity in the H dynamics. To this end, and 
also to make a direct comparison with theory, we have carried out a quantum microscopic 
calculation of the same three quantities quoted above using VMC 

The consideration of solid LiH as a quantum solid seems already justified by its Linde- 
mann ratio which is smaller than the two paradigms, '^He and ^He, but still appreciably 
larger than the common values in classical solids. From a theoretical viewpoint, it supposes 
the unavoidable introduction of at least two-body correlations to account correctly for its 
ground-state properties. We have verified using VMC that this feature also holds in solid 
LiH. 

The degree of anharmonicity in H dynamics has been established by comparing the full 
VMC calculation with another one in which the real interatomic potentials have been sub- 
stituted by harmonic approximations (Eq. (fT^ ). As commented in the previous Section, 
both VMC simulations generate identical results and then possible anharmonic effects in 
H are not observed. This situation is different from the one observed in solid ^He where 
quantum character and anharmonicity appear together. Apart from the appreciable differ- 
ence between both systems looking at their respective Lindemann ratios, a relevant feature 
that can help to understand the absence of anharmonicity of H in LiH is the significant 
difference between the interatomic potentials at short distances in both systems. Helium 
atoms interact with a hard core of Lennard- Jones type whereas the short-range interaction 
between the components of the mixture LiH is much softer (exponential type) according to 
the Born-Mayer model. 

The agreement achieved in the present work between the neutron scattering data and 
the VMC calculation is remarkably good and probably better than what could be initially 
expected from the use of semi-empirical interactions. The VMC predictions for the H kinetic 
energy and Einstein frequency coincide within error bars with the experimental measures. 
On the contrary, the H mean square displacement is about 15% larger. This points to 
probable inaccuracies of the model potentials, in particular, to core sizes smaller than the 
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real ones. In this respect, it has also to be said that these semi-empirical potentials have 
been proposed and evaluated in order to reproduce various LiH properties (e.g. lattice 
constant, bulk modulus, reststrahl frequency, etc.)i2»iS*^ in the framework of the standard 
(i.e. harmonic) lattice dynamics. Thus, the ab initio calculation of more realistic pair 
interactions could help enormously to improve a microscopic description of LiH and justify 
for the future the use of techniques beyond VMC, like the diffusion Monte Carlo or path 
integral Monte Carlo methods. 
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Li Li 7.3314 1153.80 0.0 

Li H 3.1000 187.29 0.0 

H H 5.5411 915.50 4.986 



TABLE L Parameters of the ' (r) potential by Haque and Islaniji^ 



a 


(3 










Ca/3(eVA6) 


^a/3(eVA8) 


Li 


Li 


38.48731 


17.8000 


0.2226 


0.2226 


0.0549073 


0.0216285 


Li 


H 


3.784068 


0.17675 


0.2226 


1.9000 


0.9568660 


2.0013052 


H 


H 


2.608380 


0.01422 


1.9000 


1.9000 


48.872096 


185.18284 



TABLE II: Parameters of the Vg^' (r) potential proposed by Sangster and Atwood^ estimated 
by Hussain and Sangster iil 
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FIG. 1: Neutron scattering spectra from LiH at T = 20.1(1) K: (a) TOSCA-II experimental data 
(circles) together with the estimate of the multiple scattering contribution (line); (b) experimental 
scattering law (circles) and its multiphonon component (line). 
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FIG. 2: Hydrogen-projected density of phonon states in LiH. The experimentally-determined result 
is plotted as circles, the full line represents a DDM13 lattice dynamics simulationj^^ the dashed 
line a SM7 lattice dynamics simulation/'^ and the dotted line the old neutron measurement»^ 
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FIG. 4: Variational Monte Carlo result for the H~ density profile, 'UH(r), in LiH. The solid line 
corresponds to a Gaussian with a mean squared displacement equal to the variational Monte Carlo 
value (u^) = 0.074 A^. 
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